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We have used a sample of 7 + 3 jets events collected by the DO experiment with an integrated 
luminosity of about 1 fb _1 to determine the fraction of events with double parton scattering (/dp) 
in a single pp collision at y/s — 1.96 TeV. The DP fraction and effective cross section (<r e g), a 
process-independent scale parameter related to the parton density inside the nucleon, are measured 
in three intervals of the second (ordered in pr) jet transverse momentum pi^ t2 within the range 
15 < p J y t2 < 30 GeV. In this range, /dp varies between 0.23 < /dp < 0.47, while a c g has the 
average value crjg e = 16.4 ± 0.3(stat) ± 2.3(syst) mb. 

PACS numbers: 14.20Dh, 13.85. Qk, 12.38. Qk 



I. INTRODUCTION 

Many features of high energy inelastic hadron collisions 
depend directly on the parton structure of hadrons. The 
inelastic scattering of nucleons need not to occur only 
through a single parton-parton interaction and the con- 
tribution from double parton (DP) collisions can be sig- 
nificant. A schematic view of a double parton scattering 
event in a pp interaction is shown in Fig. [TJ The rate of 
events with multiple parton scatterings depends on how 
the partons are spatially distributed within the nucleon. 



Theoretical discussions and estimations [IKE] stimulated 
measurements of DP event fractions and DP cross 
sections. The latter can be expressed as 



0"DP 



2a, 



oil' 



(1) 



where a a and up, are the cross sections of two inde- 
pendent partonic scatterings A and B. The factor m 
is equal to unity when processes A and B are indistin- 
guishable while m = 2 otherwise [1(3, EH ■ The process- 
independent scaling parameter a e g has the units of cross 
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shown in [18| that a larger fraction of DP events is ex- 
pected in the 7 + 3 jets final state as compared with the 
4-jet events. The large integrated luminosity allows us 
to select 7 + 3 jets events at high photon transverse mo- 
mentum, 60 < pi < 80 GeV (vs. p^> 16 GeV 



FIG. 1: Diagram of a double parton scattering event. 

section. Its relation to the spatial distribution of partons 
within the proton has been discussed in 0, H3, EH ■ 
The ratio ctb/^cS can be interpreted as the probability 
for partonic process B to occur provided that process 
A has already occurred. If the partons are uniformly 
distributed inside the nucleon (large a e g), ctdp will be 
rather small and, conversely, it will be large for a highly 
concentrated parton spatial density (small cr e ff). The im- 
plication and possible correlations of parton momenta 
distribution functions in ([1} are discussed in [r34l4| . 

In addition to constraining predictions from various 
models of nucleon structure and providing a better under- 
standing of non-perturbative QCD dynamics, measure- 
ments of /dp and cr e s are also needed for the accurate 
estimation of backgrounds for many rare new physics pro- 
cesses as well as for Higgs boson searches at the Tevatron 
and LHC [U,[l6|. 

To date, there have been only four dedicated measure- 
ments studying double parton scattering: by the AFS 
Collaboration in pp collisions at y/s = 63 GeV Q , by the 
UA2 Collaboration in pp collisions at yfs = 630 GeV Q , 
and twice by the CDF Collaboration in pp collisions at 
yfs = 1.8 TeV 8, 9]. The four-jet final state was used in 
the measurements to extract values of <7dp and then a e g , 
and the 7+3 jets final state was used in @ to extract /dp 
fractions and then a e g. The obtained values of a e g by 
those experiments are <7 ff w 5 mb (AFS), a e s > 8.3 mb at 
the 95% C.L. (UA2), a cS = 12.1±^° 4 7 mb (CDF, four-jet) 
and CTcff = 14.5 ± 1.7t^ mb (CDF, 7 + 3 jets). Table [Q 
summarizes all previous measurements of <7 e ff, odp, and 
/dp- 

This paper presents an analysis of hard inelastic events 
with a photon candidate (denoted below as 7) and at 
least 3 jets (referred to below as "7 + 3 jets" events) col- 
lected with the DO detector [TtJ at the Fermilab Tevatron 
Collider with yfs = 1.96 TeV and an integrated luminos- 
ity of 1.02 ± 0.06 fb -1 . In this final state, DP events 
are caused by two partonic scatterings, with 7+jets pro- 
duction in the first scattering and dijet production in the 
second. Thus, the rate of 7 + 3 jets events and their 
kinematics should be sensitive to a contribution from ad- 
ditional parton interactions. Differences in the types of 
the two final states (7+jets and dijets) and better en- 
ergy measurement of photons as compared with jets fa- 
cilitate differentiation between the two DP scatterings 
as compared with the 4-jet measurements. Also, it was 



CDF with a larger photon purity |T^|. The choice 
of a high threshold on the photon momentum provides 
(a) a clean separation between the jet produced in the 
same parton scattering from which the photon originates 
and the jets originating from additional parton scatter- 
ings and (b) a better determination of the energy scale of 
the 7+jets process. Also, in contrast to 0], the jet trans- 
verse momenta are corrected to the particle level. Other 
differences in the technique used for extracting a e s are 
described below. 

This paper is organized as follows. Section[II]briefly de- 
scribes the technique used to extract the <7 e ff parameter. 
Section IIIII provides the description of the data samples 
and selection criteria. Section IIVI describes the models 
used for signal and background events. In Section [V] we 
introduce the variables which allow us to distinguish DP 
events from other 7 + 3 jets events and determine their 
fraction. The procedure for finding the fractions of DP 
events is described in Section IVfl Section lYTll describes 
the determination of other parameters needed to calcu- 
late (T G ff. Results of the measurement are given in Section 
IVIIII with their application to selected models of parton 
density. 



TECHNIQUE FOR EXTRACTING a eB FROM 
DATA 



II. 



In the 4-jet analyses [6H8J, a g was extracted from mea- 
sured DP cross sections using Monte Carlo (MC) mod- 
eling for signal and background events and QCD predic- 
tions for the dijet cross sections. Both MC modeling and 
the QCD predictions suffer from substantial uncertainties 
leading to analogous uncertainties in <7 e ff. Another tech- 
nique for extracting cr e g was proposed in . It uses only 
quantities determined from data and thus minimizes the 
impact of theoretical assumptions. Here we follow this 
method and extract cr e ff without theoretical predictions 
of the 7+jets and dijets cross sections by comparing the 
number of 7 + 3 jets events produced in DP interactions 
in single pp collisions to the number of 7 + 3 jets events 
produced in two distinct hard interactions occurring in 
two separate pp collisions in the same beam crossing. The 
latter class of events is referred to as double interaction 
(DI) events. Assuming uncorrelated parton scatterings 
in the DP process flHa. Ill | . DP and DI events should be 
kinematically identical. This assumption is discussed in 
Appendix A. 

Measurements of dijet production with jet pr 12 — 
15 GeV [2(J in both central and forward rapidity [2l| re- 
gions indicate that the contribution from single and dou- 
ble diffraction events represents < 1% of the total dijet 
cross section. Therefore 7+jets and dijet events with jet 
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TABLE I: Summary of the results, experimental parameters, and event selections for the double parton analyses performed by 
the AFS, UA2 and CDF Collaborations. 



Experiment 


yfa (GcV) Final state p^ m (GeV) ry range 


Ceff 


C"DP, /DP 


AFS (pp), 1986 [6] 


63 


4 jets p 1 ^ > 4 


|rH<l 


~ 5 mb 


<7DP/Cr d i je t = (6 ± 1.5 ± 2.0)% 


UA2 (pp), 1991 [7] 


630 


4 jets p J T* > 15 


|r, jct | < 2 


> 8.3 mb (95% C.L.; 


a D p = 0.49 ± 0.20 nb 


CDF (pp), 1993 [8] 


1800 


4 jets p£* > 25 


|r/ jot | < 3.5 


12.1^ u 4 7 mb 


a D p = (63+^) nb, / D p = (5.4i^)% 


CDF (pp), 1997 [9] 


1800 


7 + 3 jets pj^ > 6 
p£ > 16 


|r/ jct | < 3.5 
W\ < 0.9 


14.5±1.7l2.3 mb 


/dp = (52.6 ± 2.5 ± 0.9)% 



Pt > 15 GeV are produced predominantly as a result of 
inelastic non-diffractive (hard) pp interactions. In a pp 
beam crossing with two hard collisions the probability for 
a DI event in that crossing can be expressed as 



T.IJ 



DI 



Chard Chard 



(2) 



Here a 13 and a 33 are the cross sections to produce the 
inclusive 7 + jets and dijet events, which combined give 
the 7 + 3 jets final state, and Chard is the total hard pp 
interaction cross section. The factor 2 takes into account 
that the two hard scatterings, producing a 7 + jets or di- 
jet event, can be ordered in two ways with respect to the 
two collision vertices in the DI events. The number of DI 
events, Nryi, can be obtained from Pr>i, after correction 
for the efficiencies to pass geometric and kinematic selec- 
tion criteria £di, the two- vertex event selection efficiency, 
£2vtx, and the number of beam crossings with two hard 
collisions, A^2coU : 



r 7J 



T J.I 



N m = 2- 



*2"hard <^hard 



-N2coll CDI £2v 



(3) 



Analogously to Pqi, the probability for DP events, 
Pup j in a beam crossing with one hard collision, is 



DP 



Cpp 
Chard 



Coff Chard 



(4) 



where we used Eq. ([I]). Then the number of DP events, 
Ndp, can be expressed from Prjp with a correction for 
the geometric and kinematic selection efficiency £dp , the 
single- vertex event selection efficiency ei v tx, and the num- 
ber of beam crossings with one hard collision, A^i co u: 



r 7J rr]0 



DP 



Ccff Chard 



iVicoll 6dp £h 



(5) 



The ratio of A+jp to Nr>\ allows us to obtain the ex- 
pression for a e s in the following form: 



Nbi £dp „ 

Ceff = "77 He Chard, 



(G) 



where R c = (l/2)(N lcoll /N 2coll )(e lvt Je 2vtx ). The a* 
and a 33 cross sections do not appear in this ratio and 



all the remaining efficiencies for DP and DI events enter 
only as ratios, resulting in a reduction of the impact of 
many correlated systematic uncertainties. 

Figure [5] shows the possible configurations of signal 
7 + 3 jets DP events produced in a single pp interaction 
and having one parton scattering in the final state with 
a 7 and at least one jet, superimposed with another par- 
ton scattering into a final state with at least one jet. 
We define different event topologies as follows. Events 
in which both jets from the second parton scattering are 
reconstructed, pass the selection cuts and are selected 
as the second and third jets, in order of decreasing jet 
Pt, are defined as Type I. In Type II events, the second 
jet in the dijet process is either lost due to the finite jet 
reconstruction efficiency of detector acceptance or takes 
the fourth position after the jet pr ordering. We also dis- 
tinguish Type III events, in which a jet from the second 
parton interaction becomes the leading jet of the final 
3-jets system, although they are quite rare given the pr 
range selected for the photon. 

The main background for the DP events are single par- 
ton (SP) scatterings with hard gluon bremsstrahlung in 
the initial or final state qg — > qjgg, qq — > g^gg that give 
the same 7 + 3 jets signature. They are also shown in 
Fig. [2] The fraction of DP events is determined in this 
analysis using a set of variables sensitive to the kinematic 
configurations of the two independent scatterings of par- 
ton pairs (see Sees. Mand lVT|) . 

The DI events differ from the DP events by the fact 
that the second parton scattering happens at a separate 
pp collision vertex. The DI events, with the photon and 
at least one jet from one pp collision, and at least one 
jet from another pp collision are shown in Fig. [3] with a 
similar (to DP) set of DI event types. The background to 
DI events is due to two- vertex SP events with hard 7 + 
3 jets events from one pp interaction with an additional 
soft interaction, i.e. having no reconstructed jets. The 
diagrams for these non-DI events are also shown in Fig. [3] 



III. DO DETECTOR AND DATA SAMPLES 

The DO detector is described in detail in [13] ■ Photon 
candidates are identified as isolated clusters of energy 
depositions in the uranium and liquid-argon sampling 
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DP Type II 

Jet3 
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FIG. 2: Diagrams of DP Types I, II, III and SP 7 + 3 jets 
events. For DP events, the light and bold lines correspond to 
two separate parton interactions. The dotted line represents 
unreconstructed jet. 



Dl Type I 

iet2 



Dl Type II 

Jet3 



Dl Type 




FIG. 3: Diagrams of Dl Types I, II, III and SP 7 + 3 jets 
events. For Dl events, the light and bold lines correspond 
to two separate pp interactions. The dotted line represents 
unreconstructed jet. 



calorimeter. The central calorimeter covers the pseudo- 
rapidity [22I ] range < 1.1 and two end calorimeters 
cover 1.5 < \ij\ < 4.2. The electromagnetic (EM) section 
of the calorimeter is segmented longitudinally into four 
layers and transversely into cells in pseudorapidity and 
azimuthal angle Arj x A</> = 0.1 x 0.1 (0.05 x 0.05 in the 
third layer of the EM calorimeter). The hadronic por- 
tion of the calorimeter is located behind the EM section. 
The calorimeter surrounds a tracking system consisting 
of silicon microstrip and scintillating fiber trackers, both 
located within a 2 T solenoidal magnetic field. 

The events used in this analysis should first pass trig- 
gers based on the identification of high px clusters in 
the EM calorimeter with loose shower shape require- 
ments for photons. These triggers are 100% efficient for 
p T > 35 GeV. To select photon candidates in our data 
samples, we use the following criteria [l{J. EM objects 
are reconstructed using a simple cone algorithm with a 
cone size TZ = y/ (A77) 2 + (A</>) 2 = 0.2. Regions with 
poor photon identification capability and limited p T res- 
olution, at the boundaries between calorimeter modules 
and between the central and endcap calorimeters, are 
excluded from analysis. Each photon candidate was re- 
quired to deposit more than 96% of detected energy in 
the EM section of the calorimeter and to be isolated in 
the angular region between TZ = 0.2 and TZ = 0.4 around 
the center of the cluster: (E^° t - E£° le ) / E£° TC < 0.07, 
where -E^ot is overall (EM+hadronic) tower energy in the 



matched to a reconstructed track are excluded from the 
analysis. Clusters are matched to a reconstructed track 
by computing a \ 2 function which evaluates the consis- 
tency, within uncertainties, between the reconstructed rj 
and 4> positions of the cluster and of the closest track 
extrapolated to the finely-segmented third layer of the 
EM calorimeter. The corresponding % 2 probability is 
required to be < 0.1%. We also require the energy- 
weighted EM cluster width in the finely-segmented third 
EM layer to be consistent with that expected for an elec- 
tromagnetic shower. In addition to the calorimeter iso- 
lation, we also apply a track isolation cut, requiring the 
scalar sum of track transverse momenta in a annulus of 
0.05 <TZ< 0.4 to be less than 1.5 GeV. Jets are recon- 
structed using the iterative midpoint cone algorithm [23| 
with a cone size of 0.7. Jets must satisfy quality cri- 
teria which suppress background from leptons, photons, 
and detector noise effects. To reject background from 
cosmic rays and W — ¥ tv decay, the missing transverse 
momentum in the event is required to be less than 0.7pJ. 
All pairs of objects in the event, (photon, jet) or (jet, 
jet), also are required to be separated in r\ — <f> space by 
ATZ > 0.7. 

Each event must contain at least one 7 in the rapidity 
region \y\ < 1.0 or 1.5 < \y\ < 2.5 and at least three 
jets with \y\ < 3.0. Events are selected with 7 transverse 
momentum 60 < p T < 80 GeV, leading (in p?) jet pr > 
25 GeV, while the next-to-leading (second) and third jets 
must have px > 15 GeV. The jet transverse momenta 
are corrected to the particle level. The high p T scale (i.e. 
the scale of the first parton interaction) allows a better 
separation of the first and second parton interactions in 
momentum space. 

Data events with a single pp collision vertex, which 
compose the sample of DP candidates ("lVtx" sample), 
are selected separately from events with two vertices 
which compose the sample of Dl candidates ( "2Vtx" sam- 
ple). The collision vertices in both samples are required 
to have at least three associated tracks and to be within 
60 cm of the center of the detector along the beam (z) 
axis. 

The pt spectrum for jets from dijet events falls faster 
than that for jets resulting from initial or final state radi- 
ation in the 7+ jets events, and thus DP fractions should 
depend on the jet p T 0, 110- The DP fractions 
and (T e ff are determined in three p^ t2 bins: 15-20, 20-25, 
and 25-30 GeV. The total numbers of lVtx and 2Vtx 
7 + 3 jets events remaining in each of the three p^ 2 bins 
after all selection criteria are given in Table [TTJ 

TABLE II: The numbers of selected lVtx and 2Vtx 7 + 3 jets 
events in bins of p J y t2 . 



Data 


P'r" (GeV) 




Sample 


15 - 20 20 - 25 


25-30 


lVtx 


2182 3475 


3220 


2Vtx 


2026 2792 


2309 



(77, 4>) cone of radius 1Z — 0.4 and E^° IC is EM tower en- 
ergy within a radius of 1Z = 0.2. Candidate EM clusters 
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IV. DP AND DI MODELS 

To study properties of DP and DI events and calcu- 
late their fractions in the lVtx and 2Vtx samples, re- 
spectively, we construct DP and DI models by pairing 
data events. The DP model is constructed by overlay- 
ing in a single event one event of an inclusive sample 
of 7+ >1 jet events and one event of a sample of in- 
clastic non-diffractive events selected with the minimum 
bias trig ger and a requirement of at least one jet ( "MB" 
sample) 
events. 
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Both samples contain only single-vertex 
The jet pt from the MB events is recalculated 
relative to the vertex of the 7+jet event. The resulting 
mixed events, with jets rc-ordcred in pj>, are required 
to pass the 7 + 3 jets event selections described above. 
This model of DP events, called MixDP, assumes inde- 
pendent parton scatterings, with 7 + jets and dijet final 
states, by construction. The mixing procedure is shown 
schematically in Fig. 



MixDP 




FIG. 4: Description of the mixing procedure used to prepare 
the MixDP signal sample. Two combinations of mixing 7 + 1 
jet and two jets from dijet events (a) and 7 + 2 jets and 
one jet from dijet events (b) are considered. The dotted line 
represents a jet failing the selection requirements. 

In the DI model, called MixDI, each event is con- 
structed by mixing one event of the 7+ >1 jet sample 
and one event of the >1 jet MB sample. Both events are 
exclusively selected from the two- vertices events sample. 
In the case of > 2 jets in any component of the MixDI 
mixture, the first two jets, leading in px, arc required to 
originate from the same vertex using the position along 
the beam axis of the point of closest approach to a ver- 
tex for the tracks associated to each jet and a cut on the 
minimal jet charged particle fraction, as discussed in Ap- 
pendix B. Wc consider the two- vertex 7 + jets and dijet 
events, components of the MixDI model, to better take 
into account the underlying energy, coming from the soft 
interactions of the spectator partons. The amount of this 
energy is different for single- and two- vertex events and 
causes a difference in the photon and jet identification 
efficiencies in the DP and DI events (see Section VII). 
As a background to the DI events, we consider the two- 
vertex 7 + 3 jets sample without a hard interaction at the 
second vertex (Bkg2Vtx sample), obtained by imposing 
the direct requirement that all three jets originate from 
the same vertex using the jet track information. 

The fractions of Type I (II) events in the MixDP and 
MixDI samples are the same within 1.5% for each p!^ t2 
bin and vary for both samples from 26% (73%) at 15 < 



p^ t2 < 20 GcV to (14-15)% [(84-86)%] at 25 < j^ t2 < 
30 GeV. Type III events are quite rare and their fraction 
does not exceed 1%. The MixDP and MixDI samples 
have similar kinematic (j>t and 77) distributions for the 
photon and all the jets. They differ only by the amount 
of energy coming from soft parton interactions in either 
one or two pp collisions, which may affect the photon and 
the jet selection efficiencies. 



V. DISCRIMINATING VARIABLES 

A distinctive feature of the DP events is the presence 
of two independent parton-parton scatterings within the 
same pp collision. Wc define variables sensitive to the 
kinematics of DP events, specifically to the difference be- 
tween the pt imbalance of the two object pairs in DP and 
SP 7 + 3 jets events as 0]: 



AS = A<j>{p T {l,i), p T {j,k)) 



(7) 



where the indices k (= 1,2,3) run over the jets in 
the event. Here pV(7,i) = Pt + Pt* ' an< ^ Pr{j,k) = 
p T ct ' +p T ctk , where the two object pairs, (7, jet i) and (jet 
j, jet k), are selected to give the minimal pt imbalance. 
These pairs are found by minimizing S PT , or S p i , or S<p 
defined as 



S PT 



' \Prh,i)\ 
8p T {l,i) 



\PT(j,k)\ 
Sprthk) 



1 



\Pr{l,i)\ 



\Pr{j,k)\ 



\Pt<\ 



\Pt\ 



-» j 




Pt 


+ \Pt\ 



+ 



A<j>(j,k) 



6<p(j,k) 



(8) 



(9) 



(10) 



In Eq. (flUI) A(f>(-y,i) = \ir — (f)(-f,i)\ is the supplement 
to 7T of the minimal azimuthal angle between the vectors 
Pt 7 and p T ct \ <t>(l,i)- 

The uncertainties <5p<r(7,i) in Eq. and 50(7, i) in 
Eq. (fTUj) are calculated as root-mean-square values of 
the |pr(7,*)| and A(p(j,i) distributions using the sig- 
nal MixDP sample for each of the three possible pair- 
ings. Azimuthal angles and uncertainties for jets j and k 
are defined analogously to those for the photon and jet 
i. Any of the S-variables in Eqs. ((8])— (pTjj) represents a 
significance of the pairwise p^-imbalance. On average, it 
should be higher for the SP events than for the DP events. 
Also, each S-variable effectively splits the 7 + 3 jets sys- 
tem into 7+jet and dijet pairs, based on the best pairwise 
balance. 

The two best p^-balancing pairs, which give the min- 
imum S for each of three variables in Eqs. ([5]) - (|10l) . 
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are used to calculate the corresponding AS variables, 
AS PT ,AS P ' T and AS^, according to Eq. ©. The 
AS PT , AS p ' t variables arc also used in @, , while the 
AS*0 is first introduced in this measurement. 

Figure [5] illustrates a possible orientation of the trans- 
verse momenta vectors of the photon and jets as well as 
their pt imbalances vectors, and P^, in 7 + 3 jets 
events. In SP events, the topologies with the two radia- 
tion jets emitted close to the leading jet (recoiling against 
the photon direction in </>) are preferred. The resulting 
peak at AS = it is smeared by the effects of additional 
gluon radiation and detector resolution. For a simple 
model of DP events, we have exact pairwise balance in 
Pt and thus AS will be undefined. The exact pt balance 
in the pairs can be violated due to either detector resolu- 
tion or additional gluon radiation. Both effects introduce 
an additional random contribution to the azimuthal an- 
gle between the 7+jet and the dijet px imbalance vectors, 
broadening the AS distribution (see also Fig. |H] below). 




FIG. 5: A possible orientation of photon and jets transverse 
momenta vectors in 7 + 3 jets events. Vectors and 
are the pr imbalance vectors of 7+jet and jet-jet pairs. The 
figure illustrates a general case for the production of 7 + 3 jets 
+X events. 



VI. FRACTIONS OF DP AND DI EVENTS 



A. Fractions of DP events 



where A = B1/B2 is the ratio of the background distri- 
butions, and K = (1 — /i)/(l — .{2) and C = f2l.f1 are 
the ratios of the SP and DP fractions between the DP- 
enriched and DP-depleted samples, respectively. In con- 
trast to Q, we introduce a factor A that corrects for the 
relative difference of AS shapes for the SP distributions 
in adjacent p^ t2 intervals. It is obtained using Monte 
Carlo (MC) 7 + 3 jets events generated with PYTHIA [25[ 
without multiple parton interactions and with a full sim- 
ulation of the detector response and is found to be in the 
range 0.95 — 1.3 for different bins of AS". The factor C 
is extracted using ratios of the numbers of events in data 
and MixDP samples in the adjacent bins by 



C 



^MixDP^data) j /jyMixDP/ jydataX 



(14) 



i.e. without actual knowledge of DP fractions in those 
bins. Thus, the only unknown parameter in Eq. (|13p 
is the DP fraction f\. It is obtained from a x 2 mini- 
mization of Eq. (fT3"|) using minuit [26| . The fit was per- 
formed for each pair of p 3 ^ bins (15 - 20/20 - 25 GeV 
and 20 - 25/25 - 30 GeV) and for each of AS variables 
([8"1)-(fT0]). The DP fractions in the last bin, 25 < p^ t2 < 
30 GeV, are calculated from f 2 = Cf\. The extracted 
DP fractions are shown in Fig. [51 The DP fractions, av- 
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FIG. 6: Fractions of DP events extracted with the AS 1 ^ 
AS PT , and ASy variables in the three p J T * intervals. 



To extract the fractions of DP events, we exploit 
the difference in the px spectrum of DP and radiation 
jets, mentioned in Sec. IIII1 and consider data in two 
adjacent p 1 ^ 2 intervals: DP-enriched at smaller p^, t2 and 
DP-depleted at larger p^ t2 [H, 0, H| ■ The distribution for 
each AS variable in data (D) can be expressed as a sum 
of signal (DP) and background (SP) distributions: 

Di = /iMi + (1 - (11) 

D 2 = hMi + (1 - f 2 )B 2 , (12) 

where Mj and Bi stand for the signal MixDP and back- 
ground distributions, /, is the DP fraction, (1 — /.;) is 
the SP fraction, and indices 1, 2 correspond to the DP- 
enriched and DP-depleted data sets. Multiplying (fT2"|) by 
XK and subtracting from (TTTj) we obtain: 

D 1 -XKD 2 = f 1 M 1 -XKCf 1 M 2 , (13) 



eraged over the three AS variables (with uncertainties), 
are summarized in Table Mil The location of the points 
in Fig. [6] corresponds to the mean p^, t2 for the DP model 
in a given bin. They are also shown in Table [TTT1 as (p^* 2 ) ■ 
The uncertainties are mainly caused by the statistics of 
the data and MixDP samples (used in the fitting) and 
partially by the determination of A (2 — 5)%. 

TABLE III: Fractions of DP events in the three p J ^ t2 bins. 



j eta 

Ft 


GeV 


{p>r) (GeV) 


/dp 


15- 


- 20 


17.6 


0.466 ± 0.041 


20- 


- 25 


22.3 


0.334 ±0.023 


25- 


-30 


27.3 


0.235 ±0.027 



Since each component of a MixDP signal event may 
contain two jets, where one jet may be caused by an 
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FIG. 7: Results of the two datasets fit for the AS4, vari- 
able for the combination of two p^ t2 bins 15 — 20 GeV and 
20 — 25 GeV. (a) and (b) show distributions for data (points) 
and the DP model (shaded area); (c) shows the prediction 
for DP from data (points), corrected to remove SP contribu- 
tion, and the DP model (shaded area) as a difference between 
the corresponding distributions of (a) and (b); (d) shows the 
extracted SP distributions in the two bins. The error bars 
in (a) and (b) are only statistical, while in (c) and (d) they 
represent total (statistic and systematic) uncertainty. 



additional parton interaction, the MixDP sample should 
simulate the properties of the double plus triple parton 
(TP) interactions (DP+TP), and thus the fractions in 
Table IIIII take into account a contribution from triple 
interactions as well. In this sense, the DP cross section 
calculated using Eqs. (JT|) and (|6]) is inclusive (27l. [28j. 

Figure [7] shows tests of the fit results for /1 using 
the ASj, variable for the combination of two p^ t2 bins, 
15 - 20 GeV and 20 - 25 GeV. Figure^a) show the AS 
distributions for the DP-enriched data set in data (Di) 
and the MixDP sample (Mi) weighted with its fraction 
fi . Figure EJb) shows analogous distributions for the 
DP-depleted dataset: data (D2) and the MixDP sample 
(M 2 ) weighted with its fraction / 2 . It can be concluded 
from the two distributions that the regions of small AS^ 
(< 1.5) is mostly populated by signal events with two in- 
dependent hard interactions. Figure [3c) shows the dif- 
ference between the data distributions of Figs. [3a) and 
[TJb) , corrected to remove the SP contribution by the fac- 
tor XK (the factor A corrects for the relative difference of 
the AS,p shapes and K corrects for the difference in the 
SP fractions in the two adjacent p^, t2 bins) [left side of 
Eq. (fT5|) ] and compared to the MixDP prediction [right 
side of Eq. (|13p ]. As expected, the difference is always 
positive since the fractions of DP events drop with p^ t2 . 
The DP model provides an adequate description of the 



FIG. 8: Same as in Fig. [7] but for the other combination of 
p J ^ t2 bins, 20 - 25 GeV and 25 - 30 GeV. 

data. In Fig.[7jd) we extract the SP distributions by sub- 
tracting the estimated DP contributions from the data: 
(D 1 - /!Mi)/(l - fi) for the DP-enriched data set and 
(D 2 - / 2 M 2 )/(1 - / 2 ) for the DP-depleted data sets. Fig- 
ure [8] shows the analogous test of the fit results for the 
other pair of j?° t2 bins, 20 - 25 GeV and 25 - 30 GeV. 

Predictions for SP events are obtained using pythia. 
The AS p ' t distribution for 7 + 3 jets events simulated 
with initial and final state radiation (ISR and FSR) and 
without multiple parton interactions (MPI) is shown in 
Fig. M for the interval 15 < p 1 ^ 2 < 20 GeV. Since the 
Pt imbalance of the two additional jets should compen- 
sate the pr imbalance of the "7+lcading jet" system, the 
AS„' distribution is shifted towards ir. This distribu- 
tions shows good agreement with the results for the SP 
sample shown in Fig. Ul<i). The DP 7 + 3 jets events are 
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FIG. 9: ASp^ distributions for 7 + 3 jets events simulated us- 
ing pythia with ISR/FSR but with MPI switched off (shaded 
region), as well as for 7 + 3 jets events without ISR/FSR but 
MPI switched on using Tune A-CR (triangle markers). The 
bin 15 < pfe t2 < 20 GeV is considered. 
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also simulated without ISR and FSR and using the MPI 
model corresponding to the PYTHIA parameters Tune A- 
CR [HJ. In this case, the two subleading jets may orig- 
inate only from the second parton interaction (as in DP 
events of Type I, see Fig. [2J ■ As expected, the AS P ' dis- 
tribution for these events is uniform, since the two pr 
balance vectors for the two systems, 7 + jets and dijets, 
are independent from each other. 

Another source of background to the single- vertex 7 + 
3 jets DP events is caused by double pp collisions close to 
each other along the beam direction, for which a single 
vertex is reconstructed. This was estimated separately 
and found to be negligible with a probability < 10 -3 . 



B. Fractions of DI events 

The DI fractions, /di, are extracted by fitting the 
shapes of the AS distributions of the MixDI signal and 
Bkg2Vtx background samples to that for the 2Vtx data 
using the technique described in [29|. Uncertainties are 
mainly caused by the fitting procedure and by building 
Bkg2Vtx and MixDI (in case of Type I events) mod- 
els. To estimate the uncertainty due to the Bkg2Vtx or 
MixDI models, we vary a cut on the minimal jet charged 
particle fraction (see Appendix B) from 0.5 to 0.75. The 
fitted /di in this case varies in different p 3 ^" 2 bins within 
(3 — 10)%, which is taken as the uncertainty. The final 
/di values with total uncertainties are 0.189 ± 0.029 for 
15<p^ t2 <20 GeV, 0.137 ±0.027 for 20<p£ t2 <25 GeV, 
and 0.094±0.025 for 25 <p^ t2 < 30 GeV. The relative / D i 
uncertainties grow with increasing p^ t2 . This is caused 
by a decreasing probability for a jet to originate from a 
second pp collision vertex. As a consequence, the sen- 
sitivity to DI events in the 2Vtx data sample becomes 
smaller. 

Figure [10] shows the AS^ distributions for the two- 
vertex 7 + 3 jets events selected in three p^, t2 intervals, 
15 - 20 GeV, 20 - 25 GeV and 25 - 30 GeV, for the DI 
model (MixDI) and the total sum of MixDI and Bkg2Vtx 
distributions, weighted with the DI fraction, and com- 
pared to 2Vtx data. The weighted sums of the signal 
and background samples reproduce the shapes of the data 
distributions. 



VII. DP AND DI EFFICIENCIES, R c AND a hard 

A. Ratio of photon and jet efficiencies in DP and 
DI events 

The selection efficiencies for DP and DI events enter 
Eq. (IB]) only as ratios, canceling many common correc- 
tion factors and correlated systematic uncertainties. The 
DP and DI events differ from each other by the num- 
ber of pp collision vertices (one vs. two), and therefore 
their selection efficiencies £di and £dp may differ due to 
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FIG. 10: ASr/, distributions for two-vertex 7 + 3 jets events 
in the three p£ t2 intervals: (a) 15 - 20 GeV, (b) 20 - 25 GeV 
and (c) 25 - 30 GeV. MixDI and the total sum of the MixDI 
and Bkg2Vtx distributions (shaded histograms) are weighted 
with their fractions found from the fit, compared to 2Vtx data 
(black points). The shown uncertainties are only statistical. 

different amounts of soft unclustered energy in the sin- 
gle and double pp collision events. This could lead to 
a difference in the jet reconstruction efficiencies, due to 
the different probabilities of passing the jet selection re- 
quirement pt > 6 GeV (applied during jet reconstruction) 
and different photon selection efficiencies, due to different 
amount of energy in the track and calorimeter isolation 
cones around the photon. 

To estimate these efficiencies, we use 7 + jets and di- 
jet MC events and also MixDI and MixDP data sam- 
ples. The MC events are generated with PYTHIA [25| 
and processed through a GEANT-based [3(| simulation of 
the DO detector response. In order to accurately model 
the effects of multiple proton-antiproton interactions and 
detector noise, data events from random pp crossings are 
overlaid on the MC events using data from the same time 
period as considered in the analysis. These MC events 
are then processed using the same reconstruction code as 
for the data. We also apply additional smearing to the 
reconstructed photon and jet pt so that the measure- 
ment resolutions in MC match those in data. The MC 
events arc preselected with the vertex cuts and split into 
the single- and two- vertex samples. 

The efficiencies for the photon selection criteria are 
estimated using 7 + jets MC events. We found that the 
ratio of photon efficiencies in single- vertex (e7 v ) to that 
in two- vertex samples (e^v) does not have a noticeable 
dependence on p!^ t2 and can be taken as s'{^je{ Y = 0.96+ 
0.03. The purity of 7 + jets events in the interval of 60 < 
plj, < 80 GeV in data is expected to be about 75% [l9l ]. 
and the remaining events are mostly dijet events with 
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one jet misidcntificd as photon. An analogous analysis of 
the MC dijet events gives the ratio of the efficiencies for 
jets to be misidcntificd as photons equal to 0.99 ± 0.06, 
which does not change the e{^le\ w value found with the 
signal 7 + jets sample. 

The ratio of jet efficiencies is calculated in two steps. 
First, the efficiencies are estimated with respect to a re- 
quirement to have at least three jets with p^ > 25 GeV, 
p^ t2 > 15 GeV, and p^ t3 > 15 GeV. These efficiencies 
are calculated using MC 7 4- jets and dijet events mixed 
according to the fractions of the three main MixDP and 
MixDI event types, described in Sec. IIVI The ratio of 
efficiencies for other jet selections (e.g. to get into the 
p^ t2 interval and satisfy A1Z and jet rapidity selections) 
has been calculated using MixDP and MixDI signal data 
samples. The total ratio of DP/DI jet efficiencies is 
found to be stable for all p!^ t2 bins and equal to 0.93 
with ~ 5% uncertainty. Thus, the overall ratio of photon 
and jet DP/DI selection efficiencies Edp/edi is about 0.90 
with uncertainties in the three p 3 ^ 2 bins varying within 
(5.6-6.5)%. 



B. Vertex efficiencies 

The vertex efficiency £i v tx (^2vtx) corrects for the single 
(double) collision events that are lost in the DP (DI) 
candidate sample due to the single (double) vertex cuts 
(l^vtxl < 60 cm and > 3 tracks). The ratio £i v tx/£2vtx is 
calculated from the data and found to be 1.08 ± 0.01 for 
all p 1 ^ 2 bins. The probability to miss a hard interaction 
event with at least one jet with pr > 15 GeV due to 
a non-reconstructed vertex is calculated in 7 + jets and 
minimum bias data and found to be (0.2 — 0.4)%. The 
probability to have an additional reconstructed vertex, 
passing the vertex selection requirements, is estimated 
separately using 7 + jets and dijet MC events with at 
least one reconstructed jet with px > 15 GeV and found 
to be less than 0.3%. 



C. Calculating <T h ard , N lco n and A^coii 

The numbers of expected events with one (iVi co ii) and 
two (iVicoii) PP collisions resulting in hard interactions 
are calculated from the known instantaneous luminosity 
spectrum of the collected data (Li ns t), the frequency of 
beam crossings (f CTOSS ) for the Tevatron [ljj, and the hard 
pp interaction cross section (uhard)- 

The value of Chard at y/s = 1.96 TcV is ob- 
tained in the following way. We use the inelas- 
tic cross section calculated at ^fs — 1.96 TcV, 
CT ine i(1.96 TeV) = 60.7 ± 2.4 mb [HI, found from 
avera ging the inelastic cross sections measured by the 
CDF {32[ and E811 [H Collaboration at ^/s = 1.8 TeV 



s/s = 1.96 TeV, cr SD (1.96 TcV) and cr DD (1.96 TeV), we 
use SD and DD cross sections measured at y/s = 1.8 
TeV (o- SD (1.8 TcV) = 9.46 ± 0.44 mb [H and 
cr DD (1-8 TeV) = 6.32 ± 0.03(stat) ± 1.7(syst) mb) [H[ 
and extrapolate them to -y/s = 1.96 TeV using the slow 
asymptotic behaviour predicted in [35j . We find 

<7 h axd(l-96 TeV) = CT inc i(1.96 TcV) - <r SD (1.96 TcV) 

-ct D d(1-96 TcV) = 44.76 ± 2.89 mb. (15) 

We also do analogous estimates by calculating first Chard 
at y/s = 1.8 TeV and then extrapolating it to y/s = 1.96 
TcV using EH]. This method results in er hard (1.96 TeV) = 
43.85 ± 2.63 mb which agrees well with Eq. (fT5|) . 

In each bin of the Li nst spectrum, wc calculate 
the average number of hard pp interactions (n) = 
(-^inst //cross )o-hard and then iVi co u and AT 2coll are de- 
termined from (n) using Poisson statistics. Summing 
over all L; n st bins, weighted with their fractions, we get 
Vi co n/(2V 2 coii) = 1-169 and thus i? c crhard = 56.45 ±0.88 
mb. Here we take into account that R c and Chard en- 
ter Eq. ((6]) for a e s as a product. Any increase of Chard 
leads to an increase of (n) and, as a consequence, to a 
decrease in R c , and vice versa. Specifically, while the 
found value of <7hard has a 6.5% relative uncertainty, the 
product i? c <7hard has approximately 2% uncertainty. 



VIII. RESULTS 

A. Effective cross section 

The calculation of a c g is based on Eq. ((6]) of Sec. fl] 
The numbers Ndp and A>di in each p^ 1 bin are ob- 
tained from the numbers of the lVtx and 2Vtx 7 + 3 jets 
events in Table [TTl multiplying them by /dp and /di- 
The determination of all other components of Eq. © 
are described in Sec. I VIII The resulting values of a e s 
with total uncertainties (statistical and systematic are 
summed in quadrature) arc shown in Fig. [TT] and given 
in Tabic HVl for the three p!^ t2 bins. The location of the 
points in Fig. [TT] corresponds to the mean p!^, t2 for the 
DP model in a given bin (the mean p!^ 2 values for DI 
model are the same within 0.15 GeV). These values are 
also shown in Tabic IIVI Table [V] summarizes the main 
sources of uncertainties for each p^ 12 bin. The main sys- 
tematic uncertainties are related to the determinations 
of the DI fractions (dominant uncertainty), DP fractions, 
the Edp/edi ratio, jet energy scale (JES), and i?cChardj 
giving a total systematic uncertainty of (20.5 — 32.2)%. 

TABLE IV: Effective cross section o" e ff in the three p^ t2 bins. 



and extrapolated to 1.96 TcV. To calculate single diffrac- 
tive (SD) and double diffractive (DD) cross sections at 



ict2 

Ft 


GeV 


<pO (GeV) 


a cS (mb) 


15- 


-20 


17.6 


18.2 ±3.8 


20- 


-25 


22.3 


16.3 ±3.7 


25- 


-30 


27.3 


13.9 ±4.5 
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TABLE V: Systematic (6 3y3 t), statistic (<5 s tat) and total <5totai uncertainties (in %) for o" e ff in the three p J y t2 bins. 



jetii 

F T 


Systematic uncertainty 


sources 


<5syst 


5stat 


<5total 


(GeV) 


/dp /di £dp/£di JES 


Rc Chard 


(%) 


(%) 


(%) 


15 - 20 


7.9 17.1 5.6 5.5 


2.0 


20.5 


3.1 


20.7 


20 - 25 


6.0 20.9 6.2 2.0 


2.0 


22.8 


2.5 


22.9 


25 - 30 


10.9 29.4 6.5 3.0 


2.0 


32.2 


2.7 


32.3 



The measured a c g values in the different p^ t2 bins 
agree with each other within their uncertainties, how- 
ever a slow decrease with p 1 ^ 2 can not be excluded. The 
CTcff value averaged over the three p^ 2 bins is 

all c = 16.4 ± 0.3(stat) ± 2.3(syst) mb. (16) 



£ 25r 

a" 

20 7 
15r 
10r 

5r 



D0, L = 1.0 fb 1 



20 



25 



30 



pf 2 (GeV) 



FIG. 11: Effective cross section <r c ff (mb) measured in the 
three intervals. 



B. Models of parton spatial density 

In this section we study the limits that can be obtained 
on the parameters of three phenomenological models of 
parton spatial density using the measured effective cross 
section (fT6]l . In the discussion below we follow a sim- 
ple classical approach. For a given parton spatial density 
inside the proton or antiproton p(r), one can define a 
(time- integrated) overlap 0(/3) between the parton dis- 
tributions of the colliding nucleons as a function of the 
impact parameter (3 [TTjj . The larger the overlap (i.e. 
smaller (3), the more probable it is to have at least one 
parton interaction in the colliding nucleons. The single 
hard scattering cross sections (for example, 7+ jets or di- 
jet production) should be proportional to 0(/3) and the 
cross section for the double parton scattering is propor- 
tional to the squared overlap, both integrated over all 
impact parameters f3 [H, [36| : 



Ccff 



J™0((3) 2 2tt(3 d/3 



(17) 



First, we consider the "solid sphere" model with a con- 
stant density inside the proton radius r p . In this model, 
the total hard scattering cross section can be written 



as (Thard = 47T7-2 and (T cff = Chard//- Here / is the 
geometrical enhancement factor of the DP cross sec- 
tion. It is obtained by solving Eq. (fl7|) for two overlap- 
ping spheres with a boundary conditions that the par- 
ton density p{r) = constant for r < r p and p{r) = 
for r > r p and found to be / = 2.19. The role of 
the enhancement factor can be seen better if we rewrite 
Eq. (JTJ as <tdp = /cACe/chard- The harder the singlc- 
parton interaction is the more it is biased towards the 
central hadron-hadron collision with a small impact pa- 
rameter, where we have a larger overlap of parton den- 
sities and, consequently, higher probability for a sec- 
ond parton interaction Using the measured er fF, 
for the solid sphere model we extract the proton ra- 
dius r p = 0.53 ± 0.06 fm and proton rms-radius R lms = 
0.41 ± 0.05 fm. The latter is obtained from averaging 
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r 2 as i?j? ms = J Q r 2 A-kt 2 p(r)dr = 4w J^" p(r)r 4 dr 
The results are summarized in the line "Solid Sphere" 
of Table [VlJ The Gaussian model with p(r) oc e~ r l 2a 
and exponential model with p(r) oc e~ r ^ b have been also 
tested. The relationships between the scale parameter 
(r p , a or b) and rms-radius for all the models arc given in 
Table I VII The relationships between the effective cross 
section cr c ff and parameters of the Gaussian and expo- 
nential models are taken from [38| , neglecting the terms 
that represent correlations in the transverse space. The 
scale parameters and rms-radii for both models are also 
given in Table IVII In spite of differences in the models, 
the proton rms-radii are in good agreement with each 
other, with average values varied as 0.41 — 0.47 and with 
about 12% uncertainty. On the other hand, havin g ob - 
tained rms-radius from other sources (for example, (3^1) 
and using the measured cr e fi, the size of the transverse 
correlations [3S|1 can be estimated. 



IX. SUMMARY 

We have analyzed a sample of 7 + 3 jets events col- 
lected by the DO experiment with an integrated lumi- 
nosity of about 1 fb" 1 and determined the fraction of 
events with hard double parton scattering occurring in 
a single pp collision at y/s = 1.96 TeV. These fractions 
are measured in three intervals of the second (ordered 
in pt) jet transverse momentum p^ t2 and vary from 
0.466 ± 0.041 at 15 < p^ t2 < 20 GeV to 0.235 ± 0.027 at 
25 < p 1 ? 2 < 30 GeV. 

In the same three p 3 ^ t2 intervals, we calculate an ef- 
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TABLE VI: Parameters of parton spatial density models calculated from measured o" e ff- 



Model for density p(r) 




Rims 


Parameter (fm) 


Rrms (fm) 


Solid Sphere Constant, r < r v 


A-Krl /2.2 


^3/5r p 


0.53 ± 0.06 


0.41 ± 0.05 


Gaussian e~ r2/2a2 


8na 2 




0.26 ±0.03 


0.44 ±0.05 


Exponential e~ r ^ b 


287rb 2 


sfYib 


0.14 ±0.02 


0.47 ±0.06 



fective cross section a c s, a process- independent scale 
parameter which provides information about the par- 
ton spatial density inside the proton and define the 
rate of double parton events. The measured cr c e val- 
ues agree for the three p^ t2 intervals with an average 
CT avo = i6.4±0.3(stat)±2.3(syst) mb. We note that this 
average value is in the range of those found in previous 
measurements 041] performed at different parton inter- 
action energy scales, and may indicate stable behavior of 
CTcff with respect to the considered energy scales. 

Using the measured cr e ff we have calculated scale pa- 
rameters and rms-radii of the proton for three models of 
the parton matter distribution. 
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X. APPENDIX A 

In this measurement we assume that the two parton in- 
teractions in the DP 7 ± 3 jets events can be considered 
to be independent from each other. Possible correlation 
may appear both in momentum space, since the two in- 
teractions have to share the same proton momentum, and 
at the fragmentation stage. 

In the hypothesis of two independent scatterings, the 
kinematic properties of SP dijet events should be very 
similar to those produced in the second parton interac- 
tion in the DP 7 ± 3 jets events. We compare the pt and 
rj distributions for the two cases using the pythia event 
generator, which includes momentum and flavor corre- 
lations among the partons participating in MPI. It also 
provides the possibility of choosing different MPI mod- 
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FIG. 12: Comparison of dijet events properties in SP (trian- 
gles) and in 7 + 3 jets DP events (black circles): (a) and (c) 
show comparisons of the pt and 77 distribution of the second 
(ordered in pr) jet in 7 ± 3 jets DP events with the first jet 
from the SP dijet events; (b) and (d) show comparisons of the 
Pt and 77 distribution of the third jet in 7 ± 3 jets DP events 
with the second jet from the SP dijet events. Both types of 
events are generated without ISR and FSR effects but with 
MPI Tune A-CR. 



els. In our comparison we use the PYTHIA parameters 
Tune A-CR, which is usually considered as an example 
of a model with a strong color reconnection with an ex- 
treme prediction for track multiplicities and/or average 
hadron pt [HI. As a model for the DP events, we simu- 
late 7±3 jets events using Tune A-CR but with ISR and 
FSR effects turned off and applied all selection criteria 
as described in Sec. IIII1 This configuration of the event 
generator guarantees that the two jets produced in addi- 
tion to the leading jet (and 7) in the 7 ± 3 jets final state 
arise only from additional parton interactions. The AS 
distribution for these events is shown in Fig. [9] (by tri- 
angles). The SP dijets events are also generated without 
ISR and FSR. Figure IT2Ta) compares the pr spectra of 
the first (in px) jet from the second partonic collision in 
DP events (second jet in 7±3 jets events) and the first jet 
in the SP dijet events, while Fig. [T2Tb) make analogous 
comparisons of the next (in px) jet in both event types. 
Figures [T2Tc) and [T27d) compare the 77 distributions of 
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these jets. We can see good agreement between the kine- 
matics of jets produced in the second parton interaction 
and those from the regular SP dijct events. Analogous 
comparisons were performed using Tunes A and SO with 
similar good agreement. This indicates the absence of 
visible correlations between the two DP scatterings with 
our selection criteria. 



XI. APPENDIX B 

In building signal and background DI models in 
Sec. IIV1 we take into account information about tracks 
associated with jets. We use two algorithms. In the first, 
we consider all tracks inside a jet radius (TZ = 0.7 in our 
case) and calculate the p^-weightcd position in z of all 
the tracks ("jet_z"). Here the track z position is cal- 
culated at the point of closest approach of this track to 
the beam (z) axis. For each jet in the 2Vtx data sample 
( Sec. IIIip we can estimate the distance between the jet_z 
and the pp vertex closest in z, Az(Vtx, jetj. These dis- 



tributions are shown in Fig.[l3]for each jet in the 7+3 jets 
2Vtx sample. About (95-96)% [(97-99%)] of events have 
Az(Vtx, jctj < 1.5 (2.0) cm. 

We also use an algorithm that is based on a jet charged 
particle fraction (CPF) and define a discriminant which 
measures the probability that a given jet originates 
from a particular vertex (a jet, having originated from 
a vertex, may still have tracks coming from another 
vertex). The CPF discriminant is based on the fraction 
of charged transverse energy in each jet i (in the form of 
tracks) originating from each identified vertex j in the 
event: 

EfcMtrkfSVtx,-) 



CPF(jet 4 ,Vtx,) 



EnEiMtrkTSVtx,,) 



(18) 



To confirm that a given jet originate from a vertex, we 
require Az < 2.0 and CPF > 0.5. These requirements 
being applied to two (or three) jets in the 2Vtx events 
allow to build the signal and background DI models de- 
scribed in Section IV. 
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